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ABSTRACT 

X-ray reflection spectra from photoionized accretion discs in active galaxies are presented 
for a wide range of illumination conditions. The energy, equivalent width (EW) and flux of 
the Fe Ka line are shown to depend strongly on the ratio of illuminating flux to disc flux, 
Fx/Fdisc, the photon index of the irradiating power-law, F, and the incidence angle of the 
radiation, i. When F^/F^isc^ 2 a neutral Fe Ka line is prominent for all but the largest 
values of F. At higher illuminating fluxes a He-like Fe Ka line at 6.7 keV dominates the line 
complex. With a high-energy cutoff of 100 keV, the thermal ionization instability seems to 
suppress the ionized Fe Ka line when F ^ 1.6. The Fe Ka line flux correlates with F^/F^isc^ 
but the dependence weakens as iron becomes fully ionized. The EW is roughly constant when 
Fx/Fdisc is low and a neutral line dominates, but then declines as the line progresses through 
higher ionization stages. There is a strong positive correlation between the Fe Ka EW and 
F when the line energy is at 6.7 keV, and a slight negative one when it is at 6.4 keV. This 
is a potential observational diagnostic of the ionization state of the disc. Observations of the 
broad Fe Ka line which take into account any narrow component would be able to test these 
predictions. Ionized Fe Ka lines at 6.7 keV are predicted to be common in a simple magnetic 
flare geometry. A model which includes multiple ionization gradients on the disc is postulated 
to reconcile the results with observations. 
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1 INTRODUCTION 

The discovery of the iron Ka line and Compton reflection in the 
X-ray spectra of accreting black holes was an important step in the 
unde rstanding of the central engine (M ushotzky, Done & Pounds 
1993). These features were predicted to result from the reprocess- 
ing of X-rays by an optically thick a nd cold medium (Light man & 
White 1988; [Guilbert & Rees 1988|; peorge & Fabian 199 1|; M att, 
Perola & Piro 1991), quite likely the accretion disc itself. This was 
spectacularly confirmed by the ASCA observation of a broad Fe Ka 
line in th e X-ray spectrum o f the bright Seyfert 1 galaxy MCG- 
6-30-15 ( Tanaka et al. 1995 ). The profile of this line was well fit 
by a model of line emission from relativistically moving material 
within 10 Schwarszchild radii of a supermassive black hole (F abian 
et al. 1989). Alternative explanations for such a broad line suf- 
fer from physical inconsistencies and/or fine-tuning of model pa- 
rameters (Fabian et al. 1995; Reynolds & Wilms 200C). Suddenly 
a potentially powerful probe o f accretion and bla ck hole physics 
was observationally accessible (Fabian et al. 20001. Many other ac- 
tive galactic nuclei (AGN) were subsequently observed by ASCA in 
search of a broad Fe Ka line, and, although most of the detections 



drb@ast.cam.ac.uk 



were of far less quality than the one of MCG-6-30-15, the mean 
line profile of a sample of Seyferts seemed to be broad ened in a 
similar manner (Nandra et al. 1997c ; Yaqoob et al. 2002). 



Since AGN themselves have quite variable X-ray continua, it 
was anticipated that the Fe Ka line should also change over an 
observation. The largest dataset has come from long observations 
of MCG-6-30- 15 where the line is seen to sometimes drastically 
change energy ( [wasawa et al. 1996 , 1999 ), and its flux v aries inde- 
pendently of the continuum ( Vaughan & Edelson 2001 ). This last 
propert y, also seen in other objects dWang et al. 1999T N andra et 



19997^ 



hiang et al. 2000 



JWangj^ 



Wang & Zhou 2001 



Weaver, 



Gelbord & Yaqoob 2001), is contrary to the predictions of the sim- 
ple reflection scenario, which predicts the line variations should 
track the continuum. The fact that this is not observed suggests a 
more complicate d and dynamic pattern of illumination on the disc 
( [Reynolds 2000| ) such as in the model of irradiation by m agnetic 
flares in a patchy raron^^^akev^osne^^^^ian^^^^Haardt 
& Maraschi 1 991; [Haardt, Maraschi & Ghisellini 1994[rMerlon i 
& Fabian 2( )01). These flares ir iay be rotating ( Ruszkowski 2000 ), 
outfl owing (iBeloborodov 1999|; M alzac. Beloborodov & Poutanen 
2001) and/or temporall^connecte^j^outancn^F^ 
loni & Fabian 2001), and so they will clearly have an impact on the 
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observed Fe Ka line due to the changing pattern of radiation on the 
surface of the disc. 

Many models of X-ray re flection from AGN-li ke accretion 
discs have been published (e.g. t^oss & Fabian 1993|; Zycki et a l 



Compton scattering, other emission from the gas itself, and the soft 
photons working their way outward from the dis c below) is treate d 
via the Fokker-Planck/diffusion approximation (Ross et al. 1978), 



1994; fVIagdziarz & Zdziarski 1995| ; |Ross, Fabian & Young 1999[ ). 
but they have not progressed to the same dynamic level as the mod- 
els of accretion disc coronae. The best current models compute the 
reflection spectrum from a photoionized layer on the surface of 
an accretion disc in hydrostatic equilibrium (Nayaks hin, Kazanas 
& Kall man 2000; [Ballantyne, Ross & Fabian 200 1|; R ozahska et 
al. 2001). If discs are anything like what standard theory predicts, 
then these calculations shoul d be able to make specific predictions 
about reflectioi^_s£ectr^_^ejgj_^Ja^aksl^^ & Kall- 

man 2001). 

Here, we compute such spectra from an AGN accretion disc 
over a wide range of illumination conditions, and examine the be- 
haviour of the energy, equivalent width (EW) and flux of the Fe Ka 
line. Our approach differs from that of Nayakshin & Kallman 
(2001) or [Zycki & Rozahska (2001 ) by concentrating on the ob- 
servationally accessible changes to the Fe Kq line. The parameters 
are chosen so that conditions appropriate for both flaring and qui- 
escent regions of the disc are explored. Comparison of the model 
predictions with time-averaged data (to avoid any non-equilibrium 
effects) may allow constraints to be placed on how discs are irradi- 
ated. 

The paper is structured as follows. Section ^ describes the 
model, the assumptions, and the range of parameters that are 
treated. The results of the calculations and a discussion of the evo- 
lution of the Fe Ka line profile, flux and EW are given in Section pi 
A discussion of the results is given in Section before Section g 
summarizes the primary conclusions. 



2 COMPUTATIONS 



We employ the code ofKoss & Fabian (1993 ), recently extended by 
Ballantyne et al. (200 1[), to calculate the X-ray reflection spectrum 



from the ionized layer on the surface of an AGN accretion disc. The 
deta ils of the calc ulati ons are described in the above tw o papers (see 
also Ross (1978) and Ross, Weaver & McCray 197^ ), and so only 
a brief description is given here. 

The model calculates the thermal, ionization and density struc- 
ture of the top fiveQ Thomson depths of an accretion disc that is 
irradiated by an external source of X-rays from above, and a soft 
blackbody from the disc below. The calculation is one-dimensional, 
with the heig ht of the bottom of the atmosph ere given by the disc 
equations of Merloni, Fabian & Ross (200C ) in the gas-pressure 
dominated regime and assuming a = 0.1. The flux of soft radia- 
tion entering the layer from below is assumed to result from viscous 
heating d eep in the disc and is calcul ated using the standard equa- 
tion from Shakura & Sunyaev (1973). The impinging hard X-rays 
are in the form of a power-law with photon index V (so that photon 
flux oc E~^), which strikes the surface with net flux and at an 
angle of i degrees to the normal. The incident X-rays extend from 
1 eV and terminate with a sharp cutoff at 100 keV. 

The incident radiation i s transferred analytically usi ng a one- 
stream approximation (e.g., Rybicki & Lightman 1979). Diffuse 



radiation within the layer (i.e., incident X-rays that have undergone 

^ When i^a; /iodise = 270 a total Thomson depth of 10 was used to extend 
the computational domain beyond the layer that is highly ionized by the 
illuminating radiation. 



while the Fe Kq line is treated in detail via escape a nd destruc- 
tion probabilities as suming a Lorentz line profile (see Ross 1978, 
Ross & Fabian 1993). Hydrogen and helium are assumed to be fully 



ionized everywhere within the layer, and the following ionization 
stages of the most abundant metals are treated: C V-VII, O V-IX, 
Mg IX-XIII, Si XI-XV, and Fe XVI-XXVII. The elemental abun- 



dances are given by Morrison & McCammon (1983 ). 

The major differences between this code and the one presented 
by Nayakshin et al. (2000) are in the method of radiative transfer 



(we use Fokker-Planck/diffusion, as opposed to a variable Edding- 
ton fact or scheme) and in the le vel of detail that atomic physics is 
treated ( Nayakshin et al. (200C ) employ XSTAR for the photoion- 
ization calculations and so includes many more transitions than 
our code). However, comparisons between the two programs yield 
only minor qualitative differences in the outpu t spectra, particularly 
around the Fe Ka line (Pequignot et al. 2002). 

To isolate the Fe Ka variations due to a changing illumina- 
tion pattern, the accretion disc parameters were set to represen- 
tative v alues. A typical black hole mass of 10^ Mp, ( Gebhardt et 
al. 2000) was used in the calculations, as was an accretion rate of 
0.001 of the Eddington ra te, which seems typical o f some X-ray 
bright Seyfert 1 galaxies (Morales & Fabian 2002). A radius of 
7 Schwarszchild radii was selected for where the reflection takes 
place, but we do not take into account relativistic blurring. With 
these parameters, the soft flux entering the atmosphere from below 
is -Fdisc = 1-2 X 10^^ erg cm~^ s~^, and it is assumed that there is 
no energy dissipation in the corona. Other choices for these system 
parameters would result in different Fe Ka EWs and fluxes, but 
the resulting trends with the illumination parameters should not be 
affected (Sect. 

Models were calculated for 1.5 ^ F ^ 2.1, which is the rang e 



most often observed for Seyfert 1 galaxies (Nandra et al. 1997a), 
and with the incident flux covering -Fr/-fdisc~ 0.5 - 270. To test 
for geometrical effects (which may be important for magnetic flares 
where the size of the hard X-ray emitting region is much smaller 
than the disc), for each value of F and F^c models were calculated 
for X-ray incidence angles 20 ^5 i ^5 80 degrees in steps of 20 de- 
grees. Finally, we neglect the effects of inclination along the line 
of sight, but since this does not depend on the disc or ill umination 
model, earlier results (e.g.. Matt, Fabian & Ross 1996| ) will still 
apply. 



3 RESULTS 

3.1 Fe Ka line profiles 

The properties of the Fe Ka line depend greatly o n the io nization 
state of the gas (e.g., [Matt, Fabian & Ross 1993t 1996[; R oss et 



al. 1999). When the gas is only weakly ionized or neutral (so that 
Fe I-XVI are the dominant ionization states), a strong line appears 
at 6.4 keV and the associated absorption edge is close to 7.1 keV. 
As the gas becomes more ionized, intermediate Fe states such as 
Fe XVII-XXIII are dominant. However, the Ka lines from these 
ions are preferentially suppressed due to Auger destruction (R oss, 
Fabian & Brandt 1996). It is not until the gas becomes so ionized 
that He-like iron (Fe XXV) is the dominant species that the ob- 
served Ka line shifts in energy. This ion has a large fluorescent 
yield resulting in a strong line and absorption edge at 6.7 keV and 



Fe Ka response from accretion discs 3 



~ 8.5 keV, respectively^ The line from hydrogenic Fe at 6.97 keV 
is also present in these situations, but is weaker than the helium-like 
line because of resonant scattering. Finally, it is possible for Fe to 
be fully stripped and neither a line nor an edge will be imprinted on 
the reflection spectrum. 

In variable density models it is not possible to describe the 
spectrum by assigning an ionization parameter to it. However, the 
general evolution of the Fe Ka line that was described above is 
seen in these new calculations, but the addition of a realistic density 
structure can create significant departures from this trend. Figure |l| 
shows how the Fe Ka region of the reflection spectrum can vary 
for two different values of Fx/Fg^i^c when the photon index of the 
illuminating spectrum changes. At the lower value of -Flr/Fdisc, the 
Fe Ka emission lines become more ionized (i.e., has a significant 
He-like component) with larger values of F than with smaller val- 
ues. This is possibly counter-intuitive because the harder irradiating 
spectra have more Fe-ionizing power than the softer ones. However, 
when Fj:/Fdisc is increased (the bottom panel of Fig. |l]) the situ- 
ation is reversed, and the Fe Ka line becomes highly ionized and 
weak at the lowest values of F. 

The explanation for this behaviour is illustrated in Figure ^ 
The plots show that when F = 1.5 the temperature of the photoion- 
ized gas falls abruptly from a very hot (~ Compton temperature) 
fully ionized phase to a cold, recombined phase. N ayakshin et al. 
(2000) showed that such behaviour is due to the well-known ther- 
mal ionization instabilit^^f^asesjii^ressure^alanc^^ 
McKee & Tarter 1981). There is a negligible fraction of gas with in- 
termediate ionization states, so, when F^/Fdisc is relatively small, 
only a neutral Fe Ka line is emitted with any real strength]^ At 
higher values of Fx/F^isc the incident radiation ionizes further into 
the atmosphere and the transition occurs at a non-negligible Thom- 
son depth. Few K-shell ionizing photons can traverse this layer to 
the neutral material below without being scattered outwards or ab- 
sorbed by a higher ionized species of iron. Therefore, at higher 
values of F^/Fdisc, the line weakens and disappears at the lowest 
values of P. As the photon-i ndex is increased, the maximum t em- 
perature of the gas decreases ( Guilbert, Fabian & McCray 1983 and 
Sect. |4. l| ) and the instability cannot operate. Other ionized species 
of Fe can exist near the surface and a He-like line appears in the 
spectrum. From Fig. [l] we see that one effect of this instability is to 
suppress ionized Fe Ka lines when F ^5 1.6 (but see Sect. p.l| ). 

A typical example of how the Fe Ka line evolves with increas- 
ing Fx/Fdisc is shown in Figure^. When Fx/Fdisc=0.5 or 1.0 then 
the surface of the accretion disc is not highly ionized and a neutral 
line at 6.4 keV dominates the line profile. As Fx/Fdisc is increased 
the He-like line strengthens and then weakens as the atmosphere 
is further ionized. Also, the line profile gradually broadens due to 
Compton scattering. This trend is similar to the p redictions of con- 
stant density reflection models (Ross et al. 1999). 

Finally, we will briefly discuss the impact of the incidence 
angle of the incoming radiation. When the angle of the beam is in- 
creased, the radiation encounters any specific line-of-sight Thom- 
son depth at a smaller physical depth into the atmosphere. There- 
fore, the energy is dumped into a smaller physical region, which in- 
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Figure 1. Top: The variation in the Fe Ko line profile due to changes in the 
photon index of the illuminating radiation. F increases from 1.5 at the top of 
the plot to 2.1 at the bottom, and the spectra have been vertically offset for 
clarity. Bottom: The same plot as in the other panel, but at a higher value of 
Fx/Fdisc- Here, the Fe Ka line becomes more ionized at lower F, which 
is opposite to the previous case. 



Oelgoetz & Pradhan (2001 ) showed that under certain conditions the di- 
electronic satellite lines of Fe XXV could move the energy centroid of the 
line closer to 6.6 keV. 

Thermal conduction between the two temperature zones may be impor- 
tant under some conditions, allowing enough gas at intermediate temper- 
atures to imprint spectral features on the reflection spectrum (Rozahska 



1999; U Gu&Kahn2001). 



creases the temperature of the gas (see also Nayakshin et al. 2000). 
When Fx/Fdisc is relatively small this can enhance the thermal in- 
stability, as is shown in the top panel of Figure ^. In that case, the 
He-like component in the line profile disappears when i = 80 de- 
grees. However, when Fx/Fdisc is greater (the bottom panel of 
Fig.^ the line photons are emitted at a significant Thomson depth, 
and the hot outer layers result in a weaker overall line profile. 



4 D. R. Ballantyne & R. R. Ross 




0.0001 0.001 0.01 0.1 1 

Thomson Depth 



0.1 



0.01 























r-a.o, i-60' ; 
1 1 1 



5 6 
Energy (keV) 



10 




0.0001 



0.001 



0.01 0.1 
Thomson Depth 



Figure 2. Temperature vs. Thomson Depth plots for the F = 1.5 and 1.8 
models that were shown in Fig. 0. When F = 1.5 the temperature falls 
almost non-continuously from a hot, ionized phase to a cold, recombined 
phase which is attributed to the thermal ionization instability. As Fj; /i'disc 
is increased, the transition occurs at a non-negligible Thomson depth, and 
less ionizing photons can reach neutral gas without scattering. The F = 1.8 
illuminating spectrum cannot heat the gas to as high a temperature, and so 
the instability does not occur. 



3.2 Equivalent width and line flux 

The two most important observational diagnostics of the strength 
of the Fe Kq hne are its equivalent width (EW) and the flux in 
the line. Therefore, it is important to investigate how they vary as 
the illumination of the disc is changed. To facilitate comparison 
with observations, the illuminating power-law was added to the re- 



Figure 3. The evolution of the Fe Ko line with Fx/ F^xisc for T = 2.0 and 
an incidence angle of 60 degrees. i^a;/i^disc= 0.52 for the lower spectrum 
and increases by factors of two to 270 for the upper spectrum. The spectra 
have been offset vertically for clarity. The energy evolution of the line is 
reminiscent of the one predicted by the constant density reflection models. 



flection spectrum before the measurements were made. This can 
be thought of as a reflection fraction of unity, where the reflection 
fraction measures the strength of the reflected component in an ob- 
served spectrum. In the disc-corona model of an AGN, reflection 
fractions close to one are expected. 

The EWs of the Fe Ka lines were calculated by performing 
the following integral on the total reflected-l-incident spectrum: 



EW : 



EFsjE) ~ EFcjE) 
EF^E) 



dE, 



(1) 



where EFs{E) is the spectral flux of the sum, and EFc{E) is the 
estimated spectral flux of the continuum at energy E. The contin- 
uum was estimated by fitting a straight line to the reflection spec- 
trum between 5.50 and 7.11 keV (see Figure^. The line flux was 
computed by an analogous integral over the same limits. 

The lower limit of the integration was chosen to take into ac- 
count any Compton broadening of the line profile, and the upper 
limit was set to the energy of the photoelectric absorption edge of 
neutral Fe. Using the models where Fx/Fdisc— 33, i = 20 de- 
grees and r =1.5-2.1 a comparison was made between the EWs 
calculated by this method and the ones calculated by XSPEC using 
a simple power-law plus Gaussian model. An ASCA SIS-0 response 
matrix was used to fake a 500 ks observation for each model spec- 
trum, and the fits were performed on the simulated 3-10 keV data 
with the width of the Gaussian (a) fixed at 0.25 keV. The EWs cal- 
culated by XSPEC (43-264 eV) were only ~ 50% larger than those 
obtained from the above method (26-156 eV). Larger EWs would 
result if the reflection fraction was greater than unity. 

Figure ^ shows how the flux in the Fe Kq line depends on 
the parameters of the illuminating radiation. Each panel plots the 
line flux (normalized to Fdisc) versus Fx/F^isc for the four dif- 
ferent incident angles at a specific photon index. We find that the 
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Figure 6. Plots of the Fe Kq line flux (normalized to -Fdisc) versus -FrZ-Fdisc- The separate panels differentiate between the value assumed for the photon index, 
r, of the incident power-law spectrum. As indicated, the lines in each panel distinguish between different values of the radiation incidence angle. In all cases the 
line flux correlates with the illuminating flux, but the dependence is much weaker when the line is highly ionized. 



line flux is directly proportional to the flux incident on the disc 
since the amount of energy emitted by the atmosphere must follow 
the amount that it absorbs. However, the correlation flattens as the 
line becomes highly ionized and very weak. This is noticeable for 
r ^ 1.7. An interesting feature about this plot is that it illustrates 
roughly at which -Flr/Jdisc optical depth effects become important. 
At low Fx/Fdisc the line flux increases with increasing incidence 
angle because line photons are produced closer to the surface and 
can easily escape (George & Fabian 1991). But this trend reverses 
at higher values of -Fa:/Fdisc where there is a thicker ionized skin 
and ionizing photons are easily scattered out of the layer before 
reaching any neutral material. The amount of photons "lost" to scat- 
tering increases with i and so the line flux drops as the incidence 
angle increases. The value of Fx/ F^isc at which the crossover oc- 
curs increases with F because greater fluxes are needed to produce 
a deep enough ionized layer for the softer incident spectra. 

Turning to the measurements of the equivalent width, Figure ^ 
plots the Fe Ka EW versus Fx/Fdisc- At low Fx/Fdisc a strong 
neutral Fe Ka line is present in all the reflection spectra, so the 
EWs are high. As Fx/Fdisc is increased the atmosphere becomes 
ionized, and where the ionization instability operates (F = 1.5 
& 1.6) no other line replaces the neutral line so the EW rapidly 
becomes very small. For softer illuminating spectra, a He-like line 
will replace the neutral line, and the EW can remain high, but, as 
Fx/Fdisc grows, this line will ultimately weaken and the EW will 
decrease. As seen in Fig. ^ the incidence angle crossover marks 
the onset of an optically thick ionized layer at the surface of the 
atmosphere. This point moves to higher Fx/Fdisc with larger F. 
Finally, it is worth emphasizing that for F ^ 1.8 the Fe Ka EW 
can stay roughly constant over a factor of about 10 in Fx/Fdisc- 

One other interesting correlation that can be shown with this 
dataset is illustrated in Figure ^. In this plot the models with 
Fx/Fdisc^ 2.08 are shown with thin solid lines and crosses. In 



these cases the Fe Ka line complex is dominated by the neutral 
line for most values of F, and there is only a weak dependence of 
the EW on the photon index. However, when there is a strong ion- 
ized component to the line (the curves with solid points), there is 
a significant positive correlation between the Fe Ka EW and F. 
This is a possi ble explanation for the tentativ e EW-F correlation in 
the sample of Lubinski & Zdziarski (2001). Furthermore, it may 
be possible to use these variations as an independent observational 
diagnostic to determine the ionization state of the accretion disc 
when other effects (e.g., relativistic broadening) makes it difficult 
to estimate the energy centroid of the line. 



3.3 On the self-consistency of the results 

These computations make the zeroth-order assumption that the 
changes in the irradiating flux. Fx , will not affect the accretion disc 
structure. In reality, changes in Fx are a result of the changes in disc 
structure, which will also likely affect -Fdisc- Here, we consider the 
impact of our assumption on the evolution of the Fe Ka line. 

An increase in Fx would most probably result from an in- 
crease in /, the fraction of the viscous energy dissipated in the 
corona. However, / is very likely to be a function of the accretion 
rate rh. The strength of the dependence of these three parameters 
on each other is subject to the unknown details of how energy is 
dissipated and transported in accretion discs, but the dependence 
of the disc stru cture on rh and / is well known (e.g., S vensson & 
Zdziarski 1994; [Merloni et al. 2000^ For the low accretion rate 
assumed for these calculations, the disc is gas-pressure dominated 
and the scale height of the disc H oc rn^^^{l — J)^''^". A drop 
in H will increase the density of the disc, making it more diffi- 
cult to produce ionized lines. However, this will only occur if Fx is 
strongly dependent on m and/or /. It is likely that as Fx is made 
larger by an increase in rh the disc moves from a gas-pressure dom- 
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Figure 7. Plots of the Fe Ka line equivalent width (EW) versus -FrZ-Fdisc- The sepai'ate panels differentiate between the value assumed for the photon index, 
r, of the illuminating power-law spectrum. As indicated, the lines in each panel distinguish between different values of the radiation incidence angle. For hard 
irradiating spectra, the EW will quickly decrease to very small values as the line is ionized and weakened. The decrease is postponed for softer continua where a 
strong He-like hne in the reflection spectrum will keep the EW high. 



inated phase to a radiation pressure dominated one. In this case 
H oc m(l — /), but is much larger in amplitude (by a factor of 
~ 10^) than in the gas-pressure dominated regime. Therefore, the 
density is much lower and ionized Fe Ka lines are easily produced 
over a wide range of F^/ Fdisc (e.g., F3;/Fdisc=0.48-2.88, B allan- 
tyne, Ross & Fabian 2002). 

The disc flux Jdisc oc m(l — /), and we have seen that ion- 
ized lines are common when Fx/Fdisc^ 2. Therefore, if Fx has an 
equal or greater dependence on m(l — /) then ionized lines should 
remain common. However, as commented above, there are very few 
constraints on the relations between Fx, rh, and /. Detailed mod- 
els of energy transport and dissipation in a disc-corona system are 
required to fully test the validity of these calculations. 



4 DISCUSSION 

4.1 The role of the Compton temperature 

The results presented in the previous section have shown that the 
thermal ionization instability can influence the evolution of the 
Fe Ka line in certain circumstances. This instability has been inves- 
tigated pre viously in the context of broad-line region clouds (K rolik 
et al. 1981)^n^^^^^^^^^^^^^^^^^^^^^^^^^^Jayak- 
shin et al. 2000). Here, we combine the previous work with the in- 
sight gained by our large dataset to determine, with the minimum of 
assumptions, the relative importance of the instability on the Fe Ka 
line. 



As discussed by [Krolik etal. (1981] ) and puilbert et al. (1983| ), 
the ionization instability is triggered only when the gas tempera- 
ture is very high (as seen in Figs. |l] & In a photoionized gas in 



thermal equilibrium the maximum temperature obtainable^ is the 
Compton temperature Tq: 



4fcTc / usdE - 
Ie, 



usEdE, 



(2) 



where its is the spectral energy density of the incident radiation 
field which extends from Ei to Eh- However, this equation over- 
estimates Tc when Eh > 30 keV because it does not include the 
Klein-Nishina reduction in Compton heating at higher energies. For 
radiation in the form of a single power-law ue oc _E]_£^making 



the first order Klein-Nishina correction (e.g., Ross et al. 199S) re- 
sults in 



fcTc 



(2 - V)Ei 



4(£2 



21S,e* 



(3-r) 5mec2(4-r) 



(3) 



where e = Eh/ Ei and F < 2.5. This equation can safely deter- 
mine Tc for Eh<lOQ keV. 

The dominant factors in determining Tc are e^""^ and e*""^, 
thus a larger value of T results in a lower Compton temperature. 
Alternatively, Tc rises if the irradiating spectrum is cut off at larger 
values of Ei due to the reduction in inverse Compton cooling. For 
example, with Eh = 100 keV and T = 2, fcTc = 1.3, 1.6 & 
2.1 keV at Ei = 1, 10 & 100 eV. This result and Fig. | implies 
that the lower-limit to the X-ray power-law would have to cut off in 
the hard X-rays, contrary to observations, in order to increase the 
temperature to high enough values to initiate the instability when 
r = 2. 

Other effects not described by Eq. |3| can also change Tc ■ As 



^ In reality, the maximum temperature is less than this limit due to cooling 
by bremsstrahlung radiation. 
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Figure 4. Top: The variation in the Fe Ko line profile due to changes in the 
incidence angle of the illuminating radiation, i increases from 20 degrees 
at the bottom of the plot to 80 degrees at the top, and the spectra have been 
vertically offset for clarity. Bottom: The same plot as in the other panel, but 
at higher values of F^/F^^^^ and F. The incidence angle has the strongest 
effect when Fx/F^i^^ is low, as the temperature of the gas increases with 



already noted in Sect. 3.1, Tg rises with an increase in the inci- 



dence angle (Nayakshin et al. 200C), although this can only sup- 
press ionized Fe lines when F ~ 2 if Fx / Fdisc<ti 1 where they are 
unlikely to occur for most conditions. Also, Tc can decrease if the 
disc flux increases due_^Oj_^a^j^nincrease^naccret^^ 
shin & Kallman 2001). The rise in disc flux increases the amount 



Figure 5. An example of the EW and line flux calculation. The lower spec- 
tram shows a typical computed reflection spectrum, while the upper one has 
had the illuminating power-law added to it. The continuum was estimated 
fitting a straight line between 5.5 and 7.1 keV (dotted lines). 
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Figure 8. The EW of the Fe Ka complex plotted against the photon index 
of the illuminating power-law. The thin solid lines with the crosses denote 
models with Pr/i<disc^ 2.08, while the other lines with the solid points 
(starting with the upper short-dash-long-dash line) trace out the remaining 
models with /iodise increasing downwards. When there is a significant 
ionized component to the line the correlation is opposite to the one if the 
line is predominantly neutral. 
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of Compton cooling and is, in a sense, equivalent to steepening the 
incident spectrum. 

Observations of a small sample of Seyfert 1 galaxies by Bep- 
poSAX have detected high-ener gy cutoffs in the power-law contin- 
uum between 70 and 300 keV (Matt 2000). Increasing the cutoff 



energy of the incident power-law will increase the Compton tem- 
perature at any value of F (possibly to the point of initiating the 
ionization instability; Nayakshin Sl Kazanas 2002), but it is diffi- 



cult to estimate the value of Tc for Eh > 100 keV as higher-order 
Klein-Nishina corrections need to be included. We are therefore 
placed in the interesting position where the strength of an object at 
200 keV can affect its spectral features at 6-7 keV. Models of ther- 
mal Comptonization show that the cutoff energy (~ 2kTe of the 
corona) decreases with larger F, but the exact form of the relation- 
ship depends on the geometry (e.g., Svensson 1990- If AGN with 
r=l. 8-2.1 are found to have spectral cutoffs < 200 keV then it is 
likely that these objects will exhibit ionized Fe Ka lines. However, 
it is clear that the calculations need to be extended to higher ener- 
gies to make more general predictions, and we intend to do this in 
future work. It is also necessary to obtain more measurements of 
the overall spectral shape of Seyfert 1 galaxies. 

4.2 Comparison with observations 

In principle, it is possible to compare the observed long time-scale 
variations of the broad Fe Kq line in AGN with the model pre- 
dictions. There have been many recent studies on Fe Ka variabil- 
ity based on ASCA or RXTE data (see references listed in Sect. |^ 
or, for a recent review, Reynolds 2001). However, recent XMM- 



Newton and Chandra results have shown that many AGN h ave nar- 
row Fe K q components in their spectra (|Kaspi et al. 200 1^ Y aqoob 
et al. 2001; [Reeves et al. 200 1| ; [Pounds et al. 200 1[ ). These lines are 
expected in AGN unification scenarios where reflection of X-rays 
may occur from the inner edge of a molecular torus (e.g., K rolik. 
Madau & Zycki 1994). It thus seems likely that narrow Fe Kq lines 
may be common in Seyfert galaxies. These components would have 
been unresolvable in ASCA, and so the published variations in any 
broad line would be contaminated by this feature. 

Nevertheless, it may be possibl e to compare our pred ictions to 



the observations of MCG-6-30-15. Iwasawa et al. (1999) reported 
that the broad Fe Kq line profile of this Seyfert 1 galaxy shifted 
dramatically to the red during a flare event. There was no notice- 
able narrow component {EW < 60 eV) remaining at 6.4 keV, sug- 
gesting that MCG-6-30-15 may harbour only a very weak narrow 
component. A subsequent l ong {'-^ 400 ks) o bservation of MCG-6- 
30- 1 5 by RXTE reported by Lee et al. (200C ) has resulted in a com- 
prehensive study of the Fe Kq variability. Lee et al. (200C) split 



their data into 4 different continuum flux levels and found that the 
EW of the Fe Kq line (as fitted by a Gaussian model) dropped as the 
flux increased, but that the line flux remained the same. However, 
if more complex continuum models that included reflection were 
used, the line flux decreases wi th increasing continuum flux . A re- 
analysis of this RXTE data by Vaughan & Edelson (2001) found 



that the line flux vari ed with time, but i t was not correlated with the 
continuum. As with ^ee et al. (2000| ), [Vaughan & Edelson (200 1[ ) 
found t hat the photon ind ex was correlated with the 2-10 keV flux. 
Finally, Reynolds (200C| ) used sophisticated numerical techniques 
on the RXTE data to search for a correlation between the Fe Kq line 
and the continuu m over a wider range of timescales than covered by 
Lee et al. (200C ). The search was unsuccessful which suggests that 
the line flux is uncorrelated with the continuum down to timescales 
of 0.5 ks. 
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Figure 9. The Fe Ka EW plotted against the photon index of the illumi- 
nating in'adiation, but with Fx/F^ii^^ increasing as F increases. The sohd 
points denote the evolution when Fx/F^isc= 0.52 at F = 1.5, while the 
stars show the changes when i^a;/i^disc= 4.16 at F = 1.5. In both cases. 
Fx /^disc doubles each time F increases by 0. 1 . 

Figure ^ predicts that the Fe Kq line flux should correlate 
with the flux striking the disc. Clearly, this does not agree with 
the observations. The EW was observed to drop with the contin- 
uum flux, which does agree with the predictions of Fig. Q How- 
ever, the continuum flux in MCG-6-30-15 is observed to correlate 
with the photon index, which has the opposite effect on the EW 
(Fig. Is]). To give an idea of which effect might dominate the evo- 
lution of the Fe Ka EW, we plot in Figure M how the EW varies 
with r, but tying Fx/Fdisc to the photon-index. The solid points 
shows the EW variations when Fx/Fdisc— 0.52 at F = 1.5, but 
with Fx/Fdisc doubling for an increase in F of 0.1. In this case the 
line is mainly neutral and the EW initially drops due to the increase 
in F (recall that increasing Fx/Fdisc does not change the EW by 
much when the line is dominated by a 6.4 keV line), before slowly 
increasing as the ionized component comes in. The decrease con- 
tinues the longest when i = 80 degrees. The open stars begin at 
Fx/Fdisc= 4.16 when F — 1.5 and trace the ionized evolution. 
All four panels show the same decrease-increase behaviour with 
the absolute value decreasing with incidence angle. From these 
curves and the photon-index range of MCG-6-30-15 (1.8-2.1) we 
can conclude that the EW evolution of MCG-6-30-15 is consistent 
with that of a single purely neutral line, or of a hybrid line that is ex- 
cited by radiation at a large incidence angle. However, this cannot 
explain the strange variation in line flux which may be indicative of 
more than one emission component. Indee d, both types of models 
can fit the data ( Ballantyne & Fabian 2001 ). 



4.3 Neutral and ionized Fe Ka lines 

The calculations presented here predict that ionized Fe Kq lines at 
~ 6.7 keV should be common in almost all cases if Fx/Fdiac^ 1- 
Evidence for ionized lines have been seen in the X-ray spectra 
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of so me Narrow Line Seyfert 1 Galaxies dComastri et al. 1998 
2004 [Turner, George & Nandra 1998|; [Vaushan et al. 1999|; B al 



lantyne, Iwasawa & Fabian 2001; Turner et al. 2001a | , y ) and high - 
luminosity Seyferts / low-luminosity quasars (Reeves et al. 2000 



Pounds et al. 2001 



Orr et al. 2001), both of which are thought 



to be high accretion rate objects and so might not be directly ap- 
plicable to the current computations. Most Seyfert galaxies seem 
to exhibit lines of neutral Fe at 6.4 keV (e.g. Nandra et al. 1997a), 



although this may be a result of contamination from a narrow com- 
ponent. Taking these results at face value suggests that discs are 
either weakly illuminated (i.e., Fx/Fdisc< 1) or much denser than 
previously considered. 

If discs are weakly illuminated then it is difficult to reconcile 
this with the otherwise very successful theory of coronal heating by 
magnetic flares which necessarily predicts a highly irradiated disc. 
It is possible that measurements of the centroid energy of the Fe Ka 
line can be confused by poor signal-to-noise data, calibration prob- 
lems, or gravitational redshifting effects. Searching for EW varia- 
tions with the photon index through relations shown in Fig. ^ may 
provide an alternative test for determining the ionization state of 
the Fe Ka line. 

Extensions to the basic magnetic flar e model can also recon - 
cile the models with th e data. For example, [Nayakshin et al. (2000 ), 
Nayakshin (2000) and Done & Nayakshin (2001) discuss the idea 



that a local wind will be driven by a magnetic flare event. These 
authors argue that as a result of this wind, surface material on the 
disc will be blown away revealing denser gas deeper into the disc. 
It is this gas that produces the observed reflection spectrum com- 
plete with a neutral Fe Ka line. However, the effects (or even the 
existence) of such an X-ray driven wind has not been properly mod- 
elled, so, in light of the results of this paper, the following model is 
proposed as an alternative. 

Consider the situation when separate magnetic flares are oc- 
curri ng over the accretion disc (see Figure [lo| and M erloni & Fabian 
2001). Directly underneath each illuminating event, the surface of 
the disc is highly ionized and therefore produces no Fe Kq emis- 
sion. As the angle of incidence of the radiation increases, the flux 
incident on the gas drops resulting in ionized Fe Ka lines. As the 
incidence angle further increases, the local fa; /f disc decreases un- 
til eventually neutral lines dominate. Su ch a gradient in ionization 
on the disc has also been consid ered by [Reynolds (2(300 ) for AGN 
and by [Done & Zycki (1999[), [Voung et al. (2001D and Done & 
Nayakshin (2001) to explain the low reflection fractions observed 
in low-state Galactic Black Hole Candidates. Since there is much 
more surface area on the disc that is weakly illuminated and/or illu- 
minated at grazing incidence, the majority of the line emission de- 
tected by a telescope will be at neutral energies. This effect should 
work for any value of F, as the combination of low Fx/F^iBc and 
large i will result in a neutral line (Fig. Similarly, the complete 
ionization below each flare can occur at any photon-index, but its 
assurance will depend on the physical and geometrical details of 
the event itself. 

Spectral evolution of the illumination events will also cause 
changes in the ionization state of the Fe Ka line. The recent RXTE 
observatio n of the black hole candidate GX 339-4 reported by F eng 
et al. (2001) showed that as the X-ray flux decreased from its low 
state to an "off state", its Fe Ka line evolved from neutral at 6.4 keV 
to ionized at 6.7 keV. More importantly, the photon index of the 
continuum evolved from about 1.5 to 2.2 as the flux decreased. As 
discussed earlier (Sect 4.1), when F = 1.5 the thermal ionization 
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Figure 10. A schematic picture of how ionization gradients may effect ob- 
servations of the Fe Ka line. Tliis illustration shows from a polai' viewpoint 
an accretion disc (large red to green gradient) around a black hole. On top 
of the accretion disc magnetic flares are occurring (yellow-blue regions). 
In the top half of the figure the illumination events are fairly small and lo- 
calized. Therefore, although they may highly ionize the gas directly below 
them (yellow regions), the amount of area that is producing strong ionized 
iron lines at 6.7 keV (blue regions) is relatively small compared to the area 
of the disc that is weakly illuminated. Therefore, a strong neutral line will 
be observed. This description might be appropriate for lower-luminosity ob- 
jects such as Seyferts. The bottom half of the figure shows more powerful 
and (perhaps) more numerous magnetic flares. In this case, there are large 
ai'eas of the disc that are completely ionized and so only a weak Fe Ka line 
may be detected. This could explain the X-ray Baldwin effect in quasars if 
the higher luminosity in these objects translates into a more active corona. 



served. Perhaps this is an observation of a decaying and expanding 
flare. 



4.3.1 X-ray Baldwin effect 

Observations of many Seyfert galaxies and quasars by ASCA have 
shown that as the X-ray luminosity of the source increases, the EW 
of the Fe Ka line decreases until it becomes undetectable at a 2- 
10 keV luminosity of 10^^ erg s~^ | wasawa & Taniguchi 1993| 



Nandra et al. 1997b 



Reeves & Turner 2000| ). This trend has been 



instability can operate and a neutral line is expected, but it is highly 
unlikely to work when F — 2.2, and a 6.7 keV line would be ob- 



named the X-ray Baldwin effect in analogy to the decrease in EW 
of th e C IV A1550 line in the UV continuum of QSOs (B aldwin 
1977). In terms of the constant density models, this effect could 
easily be explained as the ionization parameter of the gas increasing 
to the point where the outer layers of the disc are fully ionized 
and therefore would not imprint any line emission on the reflection 
spectrum. The ionization parameter would increase if the higher 
luminos ity in quasars was due to a higher accretion rate (M att et 
al. 1993). Newer hydrostatic disc models also can explain the X- 
ray Baldwin e ffect by ionizing a thick skin on the surfa ce of the 
accretion disk (Nayakshin 2000; Zycki & Rozahska 2001 ). 

As discussed earlier, ionization gradients along the accretion 
disc by different flare locations may be able to account for the ob- 
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served energies of the Fe Ka line. If the power and number of flares 
above the disk increases as one moves from Seyferts to quasars 
then this would increase the X-ray luminosity, and also increase 
the amount of surface area that is subject to a highly ionizing flux 
(see Fig. |l^. In such a situation only a very weak, highly ion- 
ized Fe Ka line would be emitted from the surface. Indeed, recent 
XMM-Newton observations show that ionized lines are observed in 



higher-luminosity objects (Reeves et al. 2001; Pounds et al. 2001) 



Such a scenario is not too different from that envisaged by the 
other models, but suggests a more complex geometry. There are a 
number of details to work out in order for such a model to placed 
on firm theoretical ground. Clearly, in a situation such as the one 
presented in the lower part of Fig. |l^ feedback effects between 
the reflection spectrum and the coronal activity would be vitally 
important and m ay even assist in era sing spectral features through 
Comptonization (Petrucci et al. 2001). 



5 CONCLUSIONS 

The broad Fe Ka line which (most likely) results from reflection 
from an accretion disc is the most important spectral signature in 
the X-ray study of AGN. The model calculations presented in this 
paper have shown that the behaviour and properties of the Fe line 
from an irradiated disc in equilibrium depend on how it is illumi- 
nated. Specifically, when the incident flux is larger than the disc 
flux, we predict that a He-like Fe Ka line at 6.7 keV will be a com- 
mon feature in the reflection spectrum. The thermal ionization in- 
stability is able to change this and produce a neutral line at 6.4 keV, 
but only when F ^ 1.6 (for a cutoff energy of 100 keV) and the 
incident flux is not too high. 

The Fe Ka EW will remain roughly constant as Fx/Fdisc is 
increased until the atmosphere becomes highly ionized and then 
it will decrease to small values. The rate of decline is lower for 
higher values of F. When the line complex is dominated by He-like 
Fe, the EW will strongly increase with F, but will slightly decline 
with the photon-index if it is neutral. This effect could be used as 
an observational diagnostic to determine the ionization state of the 
line, independent of the line centroid. 

Long observations of Seyferts & quasars by XMM-Newton 
would be required to test the predictions of the models. However, 
any narrow component to the line profile must be taken into ac- 
count before a comparison could be made. The results of such a 
study would be important, as it is clear that the ionization state of 
the disc depends greatly on how it is irradiated which would greatly 
constrain coronal models. 

The prediction that ionized lines are common can be recon- 
ciled with their observational rarity by employing multiple ioniza- 
tion gradients on the accretion disc, which would be a natural result 
of a magnetically active and patchy corona. The challenge is now 
extended to modellers to compute observationally testable predic- 
tions for such complex scenarios. 
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